Two-lane traffic-flow model with an exact steady-state solution 
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We propose a stochastic cellular-automaton model for two-lane traffic flow based on the misan- 
thrope process in one dimension. The misanthrope process is a stochastic process allowing for an 
exact steady-state solution; hence we have an exact flow-density diagram for two lane traffic. In 
addition, we introduce two parameters that indicate respectively driver's driving-lane preference and 
passing-lane priority. Due to the additional parameters, the model shows a deviation of the density 
ratio for driving-lane use and a biased lane-efficiency in flow. Then, a mean-field approach explicitly 
describes the asymmetric flow by the hop rates, the driving-lane preference, and the passing-lane 
priority. Meanwhile, the simulation results are in good agreement with an observational data, and 
we thus estimate these parameters. We conclude that the proposed model successfully produces 
two-lane traffic flow particularly with the driving-lane preference and the passing-lane priority. 

PACS numbers: 45.70.Vn, 05.40.-a, 89.40.-a, 47.54.-r 



I. INTRODUCTION 

Studies of traffic flow have made much progress so far 
mainly investigating highway-like single-lane traffic with 
a wide variety of theoretical and/or practical models [l|- 
7] . These models have been then extended so as to mimic 
a more realistic traffic flow, in consideration of on- /off- 
ramps, traffic light control, and multilane traffic (T), 0, 0, 
8]. In particular, two-lane traffic flow |9l-[l6j|. in which 
lane change becomes possible, is a typical extension and 
is quite meaningful in studies of traffic flow; the driver 
has a complex decision-making process in view of the 
current traffic situation when attempting lane change, 
and by which the vehicular motion is distinguished from 
the substantial motion more clearly than in the case of 
single-lane traffic. 

In most of previous works, it is a common strategy 
that one includes two general steps in the course of up- 
dating: the first step is attempt to change the lane, and 
the second is forward movement along the lane. A lane 
change is accomplished by crossing just sideways through 
"wish and decision" or when safety and incentive criteria 
are both simultaneously fulfilled; meanwhile, the forward 
movement is governed by the same dynamics as a single- 
lane model. 

In order to define a safety criterion in multilane traffic, 
they assume that the driver takes into account the dis- 
tances which the driver will have after changing the lane 
both to the vehicle ahead and behind. In some models, 
one considers the relative velocities instead of or in addi- 
tion to the distances. Anyway, the safety criterion comes 
from the first principle of traffic flow dynamics, i.e., the 
vehicle avoids a collision, and is consequently regardless 
of from the right to the left or otherwise. 

The incentive criterion reflects the second principle 
that the driver wishes to increase the velocity up to the 
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legal limit. Since lane change is possible in two-lane traf- 
fic, we can say more precisely that the driver wishes to 
avoid decreasing the velocity. This criterion is also re- 
gardless of which lane the vehicle is on. 

In addition to the above criteria, we should take legal 
restrictions into account. Particularly, a constraint on 
lane usage, which is often included in the incentive cri- 
terion, usually differs in each country as seen below [l2j : 
(a)In Germany, passing is banned in the right lane. As 
a result, they have to change from the right lane to the 
left not only when a relatively slow vehicle is ahead in 
the same lane but also when it is in the left lane. (b)In 
the United States, by contrast, passing in the right lane 
is not explicitly forbidden. We should then consider that 
the incentive criterion is symmetric to the lanes in the 
United States. (c)In Japan (where they drive on the left- 
hand side) , the left lane is assigned by law to be a driving 
lane; the right lane is, in principle, for passing. Hence, 
slow vehicles, e.g. trucks (cf. (a)), are seldom on the 
passing lane while the density of vehicles is small. 

Asymmetric rules on lane usage give rise to a differ- 
ent traffic flow in each lane, and one observes uneven 
lane use. The ratio of driving-lane use is large during 
a small number of vehicles; however it decreases quickly 
as the vehicle number increases, and falls below that of 
the passing-lane use in an intermediate density region. A 
typical inversion of the ratio for lane use is observed in 
Autobahn in Germany Also, we see an inversion in 
terms of flow as well as the density ratio in a Japanese 
expressway [l||. (As for the density ratio, one sees rather 
different behaviors each in Germany and Japan accord- 
ing to the constraint on lane usage.) Accordingly, it is a 
validation of two-lane models to reproduce an inversion 
of the density/flow ratio for lane use. 

One of the most successful models for single-lane traf- 
fic flow is the Nagel-Schreckenberg (N-S) model [l7|, a 
stochastic cellular-automaton (CA) model. Since it is ele- 
mentary as well as suitable for simulations, the N-S model 
(or its varia nts) is often incorporated as a forward-motion 
engine [l2|, [HI- Nevertheless, in comparison with the 



single-lane case, one has to pay an extra cost of renum- 
bering all the vehicles at each time step. 

In this paper, we consider a system of identical parti- 
cles hopping with nearest-neighbor interactions, focusing 
on lane-change dynamics. A complex process to try to 
change lanes is included in the hope rates which depend 
exclusively on the particle's neighboring sites. Conse- 
quently, we can specify some parameters characteristic 
to two-lane traffic flow. In Sec. II, we define the model. 
In Sec. Ill, we show the simulation results of the model. 
Analytical results corresponding to the above simulations 
are given in Sec. IV, where we consider the mean-field 
approximation as well as an exact solution. The paper 
concludes in Sec. V with some remarks. 



II. MODEL 

We propose a CA model for two-lane traffic flow based 
on processus des misanthropes [l8| : a system of identical 
particles hopping, on a lattice of finite dimension, with 
the hop rates depending not only on the occupancy of 
the departure site but also on that of the target site. 

Consider N indistinguishable particles on a one- 
dimensional periodic lattice with L sites. Each site con- 
tains at most two particles, and one site chosen randomly 
is updated at each discrete time step, i.e., the random se- 
quential updating. The particles hop from left to right 
with a given hop rates u(ni, ni+i), i.e., during the in- 
finitesimal interval dt of time, a particle hops out of 
site I occupied by ni particles into site I + 1 occupied 
by n-i+i particles with probability u(ni, ni + i)dt. (Ac- 
cordingly, the above hopping is denied with probability 
1 — u{ni,ni + i)dt.) The hop rates should take the value 
of zero except for u(l, 0), u(l. 1), u(2, 0), and u(2, 1). In 
Fig. 1(a), we illustrate a misanthrope process with hop 
rate function u(ni, This configuration is denoted 

by (m) = (0,1, 2, 1,1, 0,2, 2,1). 

We do not yet take into account which lane each parti- 
cle is on and whether they change lanes or not; only how 
many particles are contained in each road section, i.e., at 
each site of the lattice. In Fig. 1(b), we describe a con- 
figuration in the two-lane road which corresponds to Fig. 
1(a), and now the lane number is assigned to each parti- 
cle. In the following part, particles with the lane number 
(Fig. 1(b)) are referred to as vehicle in distinction from 
those in the misanthrope process (Fig. 1(a)). Provided 
that a vehicle does not cut in front of another one, we can 
specify the vehicle to hop at each time step, except for 
the configurations with the target site empty: (1,0) and 
(2,0). Accordingly, we see that the model admits two 
additional parameters, and introduce driving-lane pref- 
erence (DLP) 8 and passing-lane priority (PLP) 7. The 
DLP means the probability of a particle on the passing 
lane to attempt a lane change to the driving lane; mean- 
while, the PLP presents a priority level of the passing 
lane. In Fig. HJc), we summarize all possible motions of 
vehicles in the model and the rates to actually move. It is 
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FIG. 1. (Color online) (a)A configuration of particles in the 
misanthrope process with 9 sites and 10 particles. One parti- 
cle is now hopping from a site occupied by two to the next site 
occupied by one with rate w(2, 1). (b)A configuration of vehi- 
cles in the present model for two- lane traffic flow. We consider 
the lower cell is the driving lane and the upper cell is the pass- 
ing lane at each site. Colored cells are occupied by a vehicle; 
white cells are empty. The model coincides with the above 
misanthrope process unless one regards lanes. (c)Vehicular 
motions and the hop rates in the present model. (Any other 
motion is prohibited.) If there are two possible motions, the 
hop rate is modified with a division parameter, S or 7. 



remarkable that the present model is established without 
any conditional. 



Provided a constraint on the hop rates, 



u(2,l) = u(2,0)-tt(l,0), 



(1) 



the misanthrope process has an invariant product mea- 
sure, or in other words, an exact steady-state solution 
of factorized form [l8[ . Accordingly, as far as the mean 
values in two lanes (e.g., the total flow) are concerned, 
we can obtain exact solutions under the condition ([1}. 
The solvability condition ([1]) is derived purely through a 
mathematical argument. However, it is acceptable for the 
traffic model because generally speaking, drivers do not 
like to drive side by side, which suggests u(2, 0) > u(l, 0). 
In addition, since it cannot be comfortable to drive while 
surrounded by other vehicles, we can presume it (1,0) > 
u(2, 1). In the following part of this paper, we assume 
that the hop rates satisfy (JXJ) , giving a priority to the 
exact solution. 
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FIG. 2. (Color online) The simulation results (dot) of the proposed model for two-lane traffic with L = fOO, u(2, 0) = f., 
u(f,0) = 0.6, u(l, 1) = 0.7, u(2, 1) = 0.4, 8 = 0.02, and 7 = 0.75: (a)flow vs density, (b)velocity vs flow, (c)velocity vs 
density, (d)flow ratio for driving-lane use, and (e)density ratio for driving-lane use. The number of vehicles takes the values 
< TV < f80, changing by 5. We use 9 random configurations as the initial one. After 500 time steps, the data are collected; 
in particular, flow is measured as the average number of hops over L time steps. Also, we give analytical results corresponding 
to the above simulations: exact solutions (solid line) in (a), (b), and (c), and mean-field calculations (dashed line) in (d) and 
(e). 



III. SIMULATION RESULTS 



In Fig. 2, we show the simulation results of the present 
model: (a)flow vs density, (b)average velocity vs flow, 
(c)velocity vs density, (d)flow ratio for driving-lane use, 
and (e)density ratio for driving-lane use. Both flow and 
velocity are for two lanes. 

It is important to note that in Fig. 2, we plot the 
results with a calibration of the density The cali- 

bration method was developed to simulate particle flows 
such as vehicular traffic truly with a simple CA model. 
The cell length (i.e., site size) was changed depending 
on the vehicle density, and thus the density pcA — N/L 



used in the simulation is calibrated as 



Prw = 




where pkw means the real-world density, which is used 
in the above plots. In the calibration for single-lane 
models is given by prwv = 1 — \/l — Pca- We hence re- 
place the single-lane densities pbw & n d Pca with 2priw 
and 2pcA in order to have a calibration for two-lane mod- 
els. 

One of our goals is to compare the simulation re- 
sults with the observational data given in [l^, H(| , which 
present the characteristics of two-lane traffic flow as de- 
scribed in the Introduction. We consequently find that 
the model proposed reproduces a typical two-lane traf- 
fic flow. In particular, the inversion of the flow ratio for 
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driving-lane use (Fig. 2(d)) shows an excellent agree- 
ment with the observational data. In contrast, one does 
not see such a sharp inversion of the density (Fig. 2(e)). 
Nevertheless, we find a wide fluctuation of the ratio in a 
low density region, which also occurs in the observational 
data. It implies that the DLP and PLP surely make a 
symmetric flow unstable. 

Moreover, we consider that a strong inversion of the 
density ratio observed in real-world traffic can be at- 
tributed largely to two factors we do not take in: one is 
the existence of (relatively) slow vehicles, and the other 
is a restriction on their lane usage. In a low density re- 
gion, since slow vehicles cannot catch up with fast ones, 
slow ones hardly have an interaction with fast ones; by 
contrast, fast ones are always under the influence of slow 
ones. As a result, slow vehicles tend to get in the driving 
lane; meanwhile, fast vehicles remain in the passing lane. 



IV. ANALYTICAL RESULTS 
A. Exact solution: the flow-density plot 

As long as the hop rates satisfy the condition (fTJ), we 
can obtain an exact solution of the present model. We 
mean that from the master equation for the model, the 
steady-state probability P(oS) of finding the system in 
a configuration u> = (n±, 712, • • • , n£) is obtained in an 
explicit form and then we can calculate any expectation 
value in principle. 

In the steady state, the master equation for the mis- 
anthrope process including the present model becomes 



[w(n;_i + 1, m - l)P(. . . , + 1, m - 1, . . .) 



1=1 



(3) 



-zt(n;_i, ni)P(. . . , n h . . .) = 0, 



where < u(m, n) < 1 if m = 1, 2 and n = 0, 1; otherwise 
u(m,n) = (including the cases, m < or n < 0). 
Following [H,[25, we let 



P(w)«IJ/( 



ni ) 



(4) 



where f(ri) is called the single-site weight and one factor 
for each site of the system. Then, one immediately finds 
a solution satisfying ([3]): 



/(2) = 
"(2,1 



»(M)/(i) 2 
«(2,0)/(0) ' 
u(2,0) -u(l,0). 



(5) 



The second equation is identical to (HJ, i.e., ((T|) is a 
necessary condition. Denote by u n (uj) the number of 
sites containing n particles in a configuration lj. Since 



(Tq{uj) + <J\{io) + 0-2(10) = L and <Ji(lo) + 2o2(lj) = N, the 
first equation in leads to 



II/<",>=/<o>-/a)»(^r 



Accordingly, we have 



1 ' Z LN \u(2,0), 



where the normalization is given by 

'u{\, '\)\ CT - 



(6) 



(7) 



(8) 



The summation in @ is done over all the configurations 
lj = (n\, ri2, ■ ■ ■ , til) such that \lu\ — ni — N. One 
often refers to Z^jq as the nonequilibrium partition func- 
tion, but nevertheless the system is in a steady state but 
far from equilibrium. 

The partition function Zln is described by the Gauss 
hypergeometric function 2-F1 (a, /3, 7; z) as 



Zln — 



N 



(z£L 
\ 2 ' 



N -N+l u(l,l) 
— ,L-N + 1-A 
2 ' ' u 2,0 



. (9) 

Now, we can calculate any expectation value, in princi- 
ple, using the partition function [22|,[23|]. The flow Qln, 
defined by 



(10) 



-N 



is obtained as a function of the density p = N/L in 
the limit where L, N tends to infinity: let Q(p) — 
lim^oo Ql,pL then we obtain 



_ . u(2,0) . 
Q(p) =^r^P(2 - p) 



■-1(2.0) 



(1-P) 



(11) 



y/i-(i-4$$)p&-p)> 



1 



This gives the flow-density plot. In Fig. 2, we show 
the exact solution through the calibration ([2]). Refer to 
the Appendix „ for the detail of the calculation in this 
subsection. 



B. Mean-field approximation: lane occupancy and 
lane efficiency 



We now turn to the calculations on lane occupancy 
and lane efficiency. Here the lanes are distinguished, and 
then it seems hopeless to exactly solve the issue in the 
same manner as above. We hence use the mean-field 
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approximation in which all correlations between sites are 
neglected, i.e., we suppose 



and 



P(w) = f[n((J l ), 



(12) 



1=1 



where 7r(c) is a probability for a site of the lattice to 
be in configuration c g {e, d, p, f}; e [f] means that the 
site is empty [fully occupied], and d [p] means that one 
particle occupies the driving [passing] lane of the site. 
(Note that the distribution 7r is common in sites because 
of the periodic boundary condition.) 

Immediately one sees 7r(e) +7r(d) + 7r(p) + 7r(f) = 1 and 
p = ?r(d) + 7r(p) + 2vr(f). Then, solving © under $FZ\) 
leads to 



TT(d) 



2p 



2p + *SW(l-P+0 



5^5J-1)P(2-P) 
(13) 



J 



TT(f) 



p(2 - p) 



l)p(2-p) 



(14) 



We thus obtain 7r for a fixed p and are able to calculate 
expectation values in the mean-field approximation. 



In the approximation, the density ratio i?d of driving- 
lane use is calculated as follows: 



Rd(p) 



7T(d) + 7T(f) 




^ "(1.0) fo ~\2 



^2,0) 



(2-P) 1 



1 



(15) 



1 



<5u(l,0) 
<5u(l,0) 
""(2,0) 



«(2,0) 
(2-p) 2 + 0(2 



u(l,l)p+0(p) 2 
P) 3 . 



In Fig. 2(e), we show the graph of (|T5|) by a dashed line, 
which agrees with the simulation result. Then, from f| 15[) 
we suppose that an inversion of the density ratio does 
not occur in the present model, although the simulation 
result shows a large fluctuation of the ratio in a low den- 
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sity region. (One sees that Rd(p) is always larger than 
1/2.) Finally, expansion of in p reveals how the hop 
rates included contribute around p = 0, 2. These help us 
determine the parameter values, especially 6. 

As well, the flow ratio Rf of driving-lane use, defined 
by 



flf(p) 



u(l, 0)7r(d)7r(e) + l)7r(d)(7r(d) + tt( P )) + u{2, 0)7r(f)7r(e) + u{2, l)7r(f)7r(p) 

Q(p) 



(16) 



I 

can be calculated as a function of p. In Fig. 2(d), we expansions of R[(p) around p = 0, 2 
show the R{-Q plot using a parametric plot with p. The 
explicit expression of Rf(p) is, however, too complicated 
to read. We instead give the graph of p^l) in Fig. 3 and 



Hp) = { 



(1 + 7^(1,1) 
Su(l,0) 

_ ^ U (l,0) 
^4u(2,0) 

+ 0(2- pf. 



P + 0(p) 2 

1\ "(1, 1) S 



7~o 



2/m(2,1) 



(2-P) 



(17) 
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These also help us with parameter fitting. Working in 
the mean-field approximation, we see that the flow ratio 
has a very large fluctuation around the mean value and 
that the ratio takes similar values at the same flow but 
different densities. 
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FIG. 3. (Color online) Mean-field calculation of the density 
ratio Rf(p) for the driving-lane use. 



The results obtained may be fairly precise since in fact 
the approximation happens to show the exact solution 
pip. (It may not be so surprising but is never trivial. 
See, e.g., [25| for further discussion of mean-field theo- 
ries.) 

There are further problems to be addressed: (a)It is 
understood that the parallel update rule is better than 
the others for CA modeling of traffic flow [26J. Accord- 
ingly, we should first find an exact solution of the mis- 
anthrope process with parallel updating. (b)An open 
boundary condition should be considered from both the- 
oretical and practical viewpoints. (c)Modeling of more 
than two lane, particularly three-lane, traffic flow will be 
challenging because one is supposed to deal with a con- 
flict happening when two vehicles attempt to g et into the 
middle lane. (Then, the shuffled dynamics [27|, where ve- 
hicles are updated in random order, may be helpful.) 
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V. CONCLUSION 

In this paper, we propose a two-lane traffic-flow model 
with six rate parameters presenting free hop, overtak- 
ing, getting out of a side-by-side configuration, passing 
consecutive vehicles, the DLP, and the PLP, respectively 
denoted by u(l, 0), u(l, 1), u(2, 0), u(2, 1), S, and 7. The 
model differs substantially from previous ones in the fol- 
lowing ways: (a)In the present model, vehicles are lo- 
cated by site and lane (a so-called site-oriented descrip- 
tion) but not by the vehicle identification number (a car- 
oriented description) j24j. (b)The vehicles undertake a 
lane change either when overtaking or when returning 
to the driving lane voluntarily. Then, overtaking is cer- 
tainly done. Consequently, the DLP contributes to a 
deviation of the density ratio for driving lane use espe- 
cially in a low density region. (c)The PLP, as well as the 
DLP, explicitly gives rise to an asymmetric flow, while 
the overtaking parameter does not take the role. More 
precisely, the PLP contributes to a deviation of the flow 
ratio for driving lane use in a high density region. 

In the simulation, we optimize the parameter values in 
order to reproduce an observational data given in fl2l ]. 
and thus find that the DLP is very small and the PLP 
is relatively large. Accordingly, we think that the reason 
for changing lanes is to overtake the slower vehicle in 
front, rather than to avoid deceleration. We also take it 
that drivers have a tacit agreement that vehicles on the 
passing lane are willing to overtake those on the driving 
lane. In fact, drivers will feel unsafe if they drive side by 
side. 

Solving the master equation directly, we obtain an ex- 
act solution especially for the flow-density diagram. If 
one considers which lane a vehicle is on, the mean-field 
approach alternatively provides detailed information on 
how the parameters contribute to an asymmetric flow. 
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Appendix: Exact solution for the flow-density 
diagram 

In order to solve (O with (0| , we introduce a countert- 
erm which cancels under the sum and thus have 

/(n ; _i)/(n;) - f(ni-i)f(ni) 

= u(m_i + l,n, - l)/(n,_i + l)/(nj - l)0(m) 

(A.l) 

having cancelled common factors (a product over the 
function /(rifc) at all sites k 7^ 1—1,1), Here, f(n) is 
some auxiliary function to be determined. 

Substitution of u(m,n) = for m > 3 or n > 2 and 
f(n) = for n > 3 into (jA~T|) yields 



/(l)/(0) - /(l)/(0) = - 0)/(l)/(0), 
/(0)/(2) - /(0)/(2) = U (1,1)/(1) 2 , 
/(2)/(0) - /(2)/(0) = - u(2, 0)/(2)/(0), 
/(l)/(2) - /(l)/(2) =u(2,l)/(2)/(l). 



(A.2) 



As a result, we obtain ([5]) and 

¥M f /(0) (0 ^ WM)/(1) 2 
m <W)~ U{2 > 0) ) "(2,0)7(0) 

/(l)=(^|-u(l,0))/(l). 



(A.3) 
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Let us define A = u(l, l)/u(2,0) in the following cal- 
culation. Consider d sites are fully occupied in configu- 
ration u>, i.e., d = 02(0;) and < d < [N/2\ , then, there 
are (*^J (^Z^d) corresponding configurations. The par- 
tition function given in ([8]) is hence calculated as 



d=0 



d \N -2d 



N) A. {L -N + l) d dV 



(A.4) 



where the Pochhammer symbol (a) n = a(a + 1) • • • (a 
n — 1). Using the identity 



a\ (a + 1 



<°>- = 22 "(i)„( 



2 /n 



(A.5) 



one finds ([9]). (Note that the Gauss hypergeometric func- 
tion is defined by 2 F 1 (a, f3, T , x) = £"=0 ■) 

The flow for the system with given L and N is pre- 
sented by the partition function as 



Qln = u(u)i,u 2 )P(u)) 

\u\=N 
N 

= ^2 u ( m ' n )f( m )f( n ) 



m,n— 



Zln 



(A, 



=u(l,0) 



Zl-2,N-1 — AZ 



L-2,N-3 



2\Zl-2,N-2 + ^Zl-2,N-3 

Zln 



This is the most general form, but however, the expres- 
sion is not clear to understand. 

There is a useful identity of the Gauss hypergeometric 
function: 



(l + x) a 2 F 1 ( a ,a- 1 +l lTl x) 

/a a + 1 4x 

= 2*%,— ,7;- 



(A.7) 



2 '"(l + x) 2 /' 
by which we transform @ into a simple one, 



-N 



Z LN = [^)(l + z)- n 2 F 1 (-N,-L,L-N+l;z), 



(A.8) 



where 



A = 



(A.9) 



(A.10) 



(1 + Z)2 ' 

Then, using the following identities: 

7(7 - l)[2-Fi(o!,/3,7- l ^ x ) ~ 2Fi(a,/3,r,x)} 
— ajix 2 F\{a + + 1,7 + 1; x), 

(a - /3)(1 - x) 2 Fi(a,/3,T,x) + (7 - a) 2 Fi(a - l,P,j;x) 
+ (p-j) 2 F 1 (a,p-l, T ,x) = 0, 

(A.ll) 

x 2 F[(a,/3,T,x) = a[ 2 Fi(a + l,/3,7;x) - 2 F\{a, /3,7;x)], 

(A.12) 

we finally obtain the exact flow for arbitrary L and N, 
Qln ="(1,0)^-77 — - ( 1 - — — 



L(i-l) l-z\ N F 



u(2,0) 



N(2L-N) z / l + zF'\ 
ITl-z' ~FF 



L(L-l) 1- 



JV F. 
(A.13) 

where F = 2 Fi(—N, —L, L — N + l;z) and F' denotes 
the derivative of F with respect to z. 

In order to find the flux in the thermodynamic limit 
where L and N tend to infinity with fixed p = N/L, we 
start with the Gauss hypergeometric differential equation 
for F, 



z(l - z)F" + [1 + L-N-(1-L- N)z]F' 

-LNF = 0, 



(A.14) 



to consider g = F'/F appearing in (|A. 13|i . which is ex- 
pected to remain finite in the thermodynamic limit if 
divided by N. One sees from (|A.14[) 

z(l - z){g' + g 2 ) + [1 + L - N - (1 - L - N)z]g 

-LN = 0. 

(A.15) 

Hence we let g = g%N + g a + g-iN^ 1 + g- 2 N~ 2 H , 

finding the following quadratic equation for g\ in the iV 1 
order in the thermodynamic limit: 

pz(l - z) 9l 2 + [1 - p + (p + l)z] 9l -1 = 0. (A.16) 

Solving (|A.16|) gives 



ft = L ■ (A.17) 

Now we are ready for taking the thermodynamic limit of 
(1A~T3|) . The result is in dHJ). 
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